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ABSTRACT 

We report the discovery of a gravitationally lensed quasar identified serendipitously in the Sloan 
Digital Sky Survey (SDSS). The object, SDSS J094604.90+1 83541.8, was initially targeted for spec- 
troscopy as a luminous red galaxy, but the SDSS spectrum has the features of both iz = 0.388 galaxy 
and a z = 4.8 quasar. We have obtained additional imaging that resolves the system into two quasar 
images separated by 3.06" and a bright galaxy that is strongly blended with one of the quasar images. 
We confirm spectroscopically that the two quasar images represent a single lensed source at z = 4.8 
with a total magnification of 3.2, and we derive a model for the lensing galaxy. This is the highest 
rcdshift lensed quasar currently known. We examine the issues surrounding the selection of such an 
unusual object from existing data and briefly discuss implications for lensed quasar surveys. 
Subject headings: gravitational lensing - quasars: individual (SDSS J094604. 90+183541. 8) 



1. INTRODUCTION 

Gravitationally lensed quasars offer the opportunity to 
study quasar host galaxies at a lev el of detail not oth- 
erwise accessible (jPeng et al.ll2006D . sightlincs to study 
the physical extent of intervening absorption structures 
(Ell ison et~all[200l . and the means to constrain lens 
galaxy masses. Lensed quasar statistics are highly use- 
ful probes of cosmology, and time delays from individ- 
ual lenses provid e constraints on the Hubble constant 
(for a review see [Kochanek 2006). Several approaches 
have been adopted to construct large samples of lensed 
quasars. One method is to obtain high-resolution im- 
ages of a large number of sources to look for objects re- 
solved into multiple images on small scales. This method 
has been applied at radio wavelengths by the CLASS 
survey, from which 22 l enses were found in a survey of 
~ 16 , 000 radio sources (jMvers et al.ll2003t iBrowne et al.l 
2003). At optical wavelengths, the subarcsecond resolu- 
tion required to identify most lenses has been a chieved 
with HST snapshot surveys. iMaoz et al.1 (jl993f ) found 
five lenses in a sample of 502 quasars with z > 1; more 
recently, extensive HST follow-up of SDSS quasars at 
z > 4 yielded no le nses with > 0-1" separation i n a sam- 
ple of 162 objects ([Richards et al.ll2004bl [2006h . All of 
these searches can be considered more or less blind, in 
that little or no prior information was used to select lens 
candidates. 

An alternative is to select likely candidates for lensed 
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quasars directly from imaging data. This is the ap- 
proach taken by the SPSS Quasar Lens Search (SQLS: 
lOeuri et al.ll2006l 120081 llnada et al.l [2008 and references 
therein. iKayo et al.l l2009h . which has identified 40 new 
gravitational lenses through follow-up of quasars within 
the SDSS. Candidates are selected either by the presence 
of a nearby object with similar colors as the quasar, or 
by the quasar having image being resolved, a hint that 
it may represent the blend of multiple images at small 
separations. 

Only a handful of z > 2.5 quasars are currently known 
to be gravitationally lensed. Both blind and targeted 
surveys for lensed high redshift quasars are limited by 
the extreme rarity of the source po pulation. That no 
lenses were found in the survey of iRichards et al.l - 
which included more than half of the z > 4 quasars 
known at the time — is unsurprising given lensing rates 
of < 1% found by surveys of quasars a t lower red- 
shifts (|Inada et al.ll2008t IBrowne et al.ll200l); though the 



lensin g probability increases with redshift (| Turner et al.1 
H98l and magnification bias is expected to further en- 
hance the observed lensing rate at high redshift in flux- 
limited surv eys (for discussions o f thi s effect on SPSS 
quasars, see iWvithe fc Loebl 12001 and iComerford eTaLI 
2002). Targeted surveys such as the SQLS further de- 
pend on the initial selection method employed to iden- 
tify high-redshift quasar candidates. In the case of the 
SPSS, quasar targets at z > 2.5 are required to have a 
point-like morphology, introducing a strong bias against 
small-separation lensed systems. 

We present the discovery of a z = 4.8 lensed quasar 
drawn from the SPSS. This is the highest redshift 
lensed quasar currently known, and only the third 
at z > 4 after BRI 0952-011 5 at z = 4.426 (1" 
separation, iMcMahon et al.l I1992D a nd PSS 2322+1944 
at z = 4.12 (1.4 9" separation, iCarilli et all 120021: 
IRiechers et aLll2008h . SPSS J094604.90+183541.8 (here- 
after J0946+1835) was initially targeted for spectroscopy 
as a luminous red galaxy and identified by strong quasar 
emission lines apparent in the SPSS spectrum. 
This paper is organized as follows: in Section [2] we 
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Fig. 1.— SDSS color image of J0946+1835. The image is cen- 
tered on the photometric object targeted as a galaxy; the strong 
blending from additional, redder objects is apparent. 



TABLE 1 
SDSS Photometry 



Object 



A > 23.4 
B+G 20.90 ± 0.06 
GALFIT results: 
A' > 23.4 
B' > 23.4 
G' 20.89 ±0.21 



21.34 ±0.14 
18.87 ±0.01 

21.42 ± 0.04 
20.46 ± 0.05 
18.95 ± 0.06 



19.52 ±0.12 
17.98 ±0.01 

19.49 ±0.02 
19.24 ±0.04 
18.30 ±0.04 



19.46 ±0.11 
17.66 ± 0.02 

19.53 ±0.05 
19.30 ±0.17 
17.80 ±0.08 



Note. - The SDSS positions (J2000) are 09:46:04.90 
+ 18:35:41.8 for A and 09:46:04.79 ±18:35:39.7 for B+G. For ob- 
ject A the values are from PSF magnitudes, while for object B+G 
the values are model magnitudes; the S DSS photometry is given 
in asinh magnitudes (Lupton ct al. 1993). Deblended photom etry 
from GALFIT is given in the last three rows and described in i|2.2l 
magnitudes are on the AB system. 



show how the object was serendipitously identified from 
SDSS data and provide results from extensive photomet- 
ric and spectroscopic follow-up that confirm the lensing 
hypothesis. We then provide results from modeling the 
lens galaxy (§ [3]), briefly discuss the nature of the dis- 
covery of this object (§H]), and conclude with prospects 
for future searches for high-redshift quasar lenses (§[5]). 
Unless otherwise noted, all magnitudes are on the AB 
system and corrected for Ga lactic reddening using the 
maps of lSchlegel et al.l ((19981 ). We adopt a standard cos- 
mology with parameters Ho = 70 km s _1 Mpc^ 1 , Q m = 
0.3,fi A = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 

J0946+1835 was targeted for spectroscopy by the 
SDSS as a Luminous Red Galaxy (LRG) according to the 
qri co lor selection technique outlined inlEisenstci n et al.l 
pooil) . As we will discuss in § 2J this selection occurred 
because the foreground galaxy is significantly brighter 
than the background quasar images. In this section we 
present various observations of this system, beginning 
with the SDSS imaging and spectroscopy from which it 
was initially identified, followed by higher resolution and 
signal-to-noise ratio imaging, and finally spectroscopic 
observations of the multiple components comprising the 
system. All of these data provide strong evidence in sup- 
port of the lensing hypothesis. 
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Fig. 2.— SDSS spectrum of J0946+1835, smoothed with a 5 
pixel boxcar function. The Lyc* and N V emission lines unambigu- 
ously show the signature of a z = 4.8 quasar. However, the level 
of flux blueward of Lyo (and in particular blueward of the Lyman 
limit) indicates excess light from another object; an expansion of 
the region around 5500A shows that this object is a z = 0.388 
galaxy. 



2.1. SDSS Imaging & Spectroscopy 

The SDSS is a photometric and spectro scopic sur- 
vey c overing nearly one-quarter of the sky (lYork et al 
2000). Images with the multi-CCD camera (jGunn et al 
19981 ) are obtained in five optical bands (ugriz) and re- 



duced accord i ng to a photometric pipeline described in 
iLupton et al.l (|2001l ). Details of the photome t ric ca libra- 
tion procedure are given in iFukugita etUI JT996) and 
(12002D. and the astrometric calibration is 
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outlined in iPier et al.l ()2003f ). Photometric objects are 
selected for spectroscopy by y arious programs as de- 
tailed in iStoughton et a l. (2002) , including, amo ng oth- 
ers, ~ 10 5 quasar targets ([Richards et al.l l2~002af ) and a 
similar number of Lum inous Red Galaxy (LRG) targets 
(jEisenstein et al.ll200lD . 

J0946+1835 was imaged during the course of the SDSS 
main survey on 2005 Mar 10 with seeing of 1.2", 1.2", 
and 1.5" in the r, i, and z bands, respectively. The SDSS 
photometric pipeline identified two objects in this field: a 
stellar object, and a much brighter galaxy located ~ 2.4" 
to the SW. As we will show, this galaxy is actually an 
unresolved blend of a low-redshift galaxy and a high- 
redshift quasar. Although the galaxy is contaminated 
with quasar light, it was selected as an LRG spectro- 
scopic target based on its gri colors. De tails of the SPSS 
photo metry from Data Release 7 (DR7, [Abazajia n et al.l 
2009]) for the two objects are given in Table [TJ 

SDSS spectra are acquired with a multi-fiber spectro- 
graph, with each fiber subtending a 3" diameter region 
on the sky. The spectra span 3800-9200A at a resolution 
of ~ 1800; the wavele ngth and flux ca l ibrati on of the 
spectra are described in Stoughton ct al. (2002). A spec- 
trum was obtained on 2006 Jan 24 with 1.3" seeing. The 
SDSS fiber was centered on the position derived for the 
B+G object. The spectrum clearly shows the emission 
line features of a z = 4.8 quasar (Figure [5]). The spec- 
trum also shows the Ca H+K absorption features of a 
z = 0.388 galaxy. The automated spectroscopic pipeline 



A Lensed Quasar at z = 4.8 



3 



5" 





B 




2 PSFs + 1 G sub 

Fig. 3.— ARC SPIcam imaging of J0946+1835 with 0'.'9 see- 
ing. North is up and east is left. The upper left panel shows a 
cutout image from the i band. Three components have been fit 
with GALFIT; the upper right panel shows the residual image af- 
ter the model galaxy has been removed, while the lower left panel 
shows the residual image after the two model PSFs have been re- 
moved. The lower right panel shows the residual image with all 
model components subtracted, with lines whose intersections in- 
dicate the positions of each of the components. A faint galaxy 
[i = 22.1) remains visible to the NW, as well as a faint excess to 
the east. 

spectrold (|Stoughton et al.l I20Q2D classified the object 
as a quasar with low confidence (56%) and assigned a 
redshift z = 4.806. 

The complex nature of the spectrum became appar- 
ent from manual inspectio n during construction of the 
SDSS DR7 quasar catalog ([Schneider et all 120101 ) . One 
author (DPS) identified 1083 spectra from DR6 and DR7 
with quasar classifications but with questionable red- 
shifts. From among those objects another author (PBH) 
identified six quasar/galaxy spectral blends whose SDSS 
images were consistent with a lensing hypothesis. Of 
those six, one had already been been studied in the SQLS 
and found to be a quasar-galaxy pair, one was deter- 
mined to be a quasar next to a z — 0.0792 galaxy, three 
are being followed up further, and the sixth is the object 
reported here. 

2.2. ARC Imaging 

J0946+1835 was observed with the SPIcam imager on 
the ARC 3.5m telescope at Apache Point Observatory 
on 2009 Nov 6. The 2048 2 backside-illuminated SITe 
CCD was binned by a factor of two along each axis dur- 
ing readout, resulting in a pixel scale of 0.28" over the 
4.8' field of view. The seeing was 0.9" and conditions 
were photometric. Two dithered exposures of 240s each 
were obtained in the i-band and two dithered exposures 
of 300s each were obtained in the z-band. The images 
were reduced using standard IRAF0 routines called from 

8 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
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Fig. 4. — Contour plots of the ARC i-band imaging. Left panel: 
Total flux contours. Right panel: Contours of flux for the fitted 
components, including the two PSFs (black lines) and the model 
galaxy (gray lines). All contours are drawn at (19.6, 20.2, 21.4, 
22.6) mag arcsec -2 , and positions are relative to A. 

Pyrafl scripts. After applying a bias correction, super- 
sky images for both filters were constructed by combining 
science images from other fields obtained just before the 
observations of J0946+1835; a total of seven i-band and 
23 z-band images were used to construct the flats. Both 
astrometric registration and photometric calibration of 
the images were accomplished using bright SDSS stars 
within the field. 

We derived positions and fluxe s for the various c ompo- 
nents using the GALFIT code (jPeng et all 120021 For 
both the i and z images, approximate positions and 
fluxes obtained from examining the images with IRAF 
were used as initial conditions for the fitting. An empir- 
ical PSF was constructed from a set of 10 bright stars 
within the field. The fitting was then performed using 
two PSFs to represent the quasar images and a Sersic 
model for the galaxy. Results of the fitting are shown in 
Figure [3] and Tabled For the i-band image, the separa- 
tion between the quasar images is 3.06" and the flux ratio 
is fs/fA = 0.7. Isophotal contours for the model compo- 
nents derived from the i-band imaging are displayed in 
Figure SI showing the strong blending between B and G. 
Table [3] shows parameters for the Sersic profile derived 
by GALFIT for the lens galaxy. 

A faint galaxy is also seen ~ 5" to the northwest (la- 
beled G2 in Figure |3j . This galaxy is not detected by 
the SDSS pipeline, but we were able to extract aper- 
ture photometry in the g, r, and i bands, and obtain a 
flux of r — 22.0 ± 0.1 and colors g — r = 1.3 ± 0.3 and 
r — i = 0.3 ± 0.2. Photometry from the ARC imaging 
yields i = 22.1 ± 0.1 and i - z = 0.3 ± 0.2. These col- 
ors are inconsistent with an object at the redshift of the 
quasar, indicating that it is not an additional lensed im- 
age. However, this faint galaxy may be associated with 
the bright lensing galaxy. There also appears to be a flux 
excess (considerably fainter than galaxy G2) to the east 
of quasar image A in the SPIcam i band data. 

2.3. UH88 Imaging 

Additional imaging data of J0946+1835 were obtained 
on 2010 Jan 16 using the University of Hawaii 2.2m tele- 
scope (UH88) with the Tektronix 2048 2 CCD camera, 

for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

9 Pyraf is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 
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TABLE 2 

Relative Astrometry and Photometry 











UH88 




ARC 


Component 


AX (") 


AY (") 


V 


R 


1 


i 


z 


A 


0.000 ± 0.003 


0.000 ± 0.003 




21.05 ±0.02 


19.90 ±0.01 


19.76 ± 0.01 


19.71 ±0.02 


B 


1.756 ±0.005 


-2.503 ± 0.005 




21.74 ±0.07 


20.00 ±0.04 


20.12 ±0.03 


19.82 ±0.05 


G 


1.453 ± 0.005 


-1.942 ±0.008 


19.37 ±0.06 


18.75 ±0.02 


18.06 ±0.02 


18.34 ±0.01 


17.96 ± 0.02 



Note. — The reference position for component A is 09:46:04.90 ±18:35:41.85 (J2000), with a systematic error in mapping to SDSS 
astrometry of 0'.'06 in each coordinate. All positions are derived from the ARC i-band imaging; positive X and Y represent west and north, 
respectively. Magnitudes are based on PSF models for A and B and Sersic model for G, as calculated by GALFIT. 



TABLE 3 
Parameters of the galaxy profile 



V R I i z 



Re(") 2.59 ±0.26 1.56 ±0.07 1.89 ±0.11 1.24 ± 0.03 1.33 ± 0.05 

e 0.29 ±0.05 0.16 ±0.02 0.13 ±0.02 0.13 ±0.01 0.13 ±0.02 

9 e {°) 28.8 ±6.5 14.2 ±5.4 26.7 ±7.0 16.6 ± 5.0 24.9 ± 8.6 

Note. — All measurements taken from GALFIT results based on fitting a Sersic profile with n = 4. The parameters are the effective 
radius (Re), ellipticity (e = 1 — b/a), and position angle of the ellipticity (9 e , measured east of north). 



TABLE 4 
Multi-epoch i-band Photometry 



Source 


Epoch 


A 


B 


SDSS 


2005.19 


19.49 


19.24 


ARC 


2009.85 


19.76 


20.12 


UH88 a 


2010.04 


19.82 


19.92 



Note. — Photometry for the two quasar images over three 
observational epochs. 



a UH88 i-band magnitudes are derived from the /-band imaging 
after computing a color difference from the observed SED of I — i = 
0.08. 

with a pixel scale of 0.22". Conditions were clear and 
photometric with 0.8-1.0" seeing. Images were obtained 
in the VRI bands; each band received four exposures 
of 120s each. Images were reduced using standard IRAF 
routines and photometric calibration was provided by the 
SA92 field. Magnitudes were derived on the Vega system 
and converted to AB using Vab = VVega + 0.00, Rab = 
i?vc ga + 0.20, and 7 A b = ^Vcga + 0.45. 

The UH88 images were also fit by three components us- 
ing GALFIT using a similar process to that of the ARC 
imaging. Results of fitting the UH88 images are given in 
Tables [5] and [3] Overall, there is a high level of consis- 
tency between the positions derived by GALFIT in all 
of the imaging bands, and also between the UH88 and 
ARC imaging. 

The galaxy G2 is detected in V, R, and /. As with the 
SPIcam i-band image, excess flux appears to the east of 
quasar image A in the UH88 / band image, indicating a 
real detection. There is a ~ 1.4er excess at this location in 
the V band; this object is also likely to be a low-redshift 
galaxy. 

2.4. Imaging Summary 

We have presented results from three epochs of imag- 
ing for this system, summarized in Table There are 
some apparent differences when comparing the observa- 
tions. We first note that the quoted photometric errors 
(Tables [1] and [2]) are purely statistical errors from GAL- 



FIT and do not reflect the uncertainties arising from the 
fitting itself and from the photometric calibration. The 
fitting uncertainties are particularly important for the 
SDSS imaging, which has the poorest resolution and the 
largest pixel size, hampering robust image decomposi- 
tion. We have derived deblended photometry from the 
SDSS images by using GALFIT with some constraints to 
improve the fitting, namely, we fix the positions of the 
three components to the ARC-derived positions, and fix 
th e g alaxy fluxes to be within 0.1 mag of the ARC val- 
ue£j. The results of this fitting are shown in Table [TJ 
We also experimented with unconstrained fitting of the 
SDSS images and had poor results. In particular, offset- 
ting the initial positions by just a few pixels generally 
results in a fit with only two components, component 
B effectively being lost. This highlights the difficulty 
of fully characterizing a system such as this from SDSS 
imaging alone. 

The colors of component B are slightly redder than 
those of component A in both the ARC and UH88 imag- 
ing, but not in the SDSS imaging. The ARC colors are 
(i ~ z) A = 0.05 ± 0.02 and (i - z) B = 0.30 ± 0.06; on 
the other hand, the SDSS colors are i ~ z w for both 
components. The redder color of component B may be 
due to greater extinction from the lens galaxy; however, 
we argue against differential reddening of the two images 
in § 12 51 Some of the observed color differences between 
the two components are due to emission line differences 
of the quasar spectra, discussed in § 12.51 The remainder 
are likely due to the difficulty of decomposing the image 
properly, and are indicative of systematic errors in the 
photometry greater than the quoted statistical errors. 

Component B is significantly brighter in the SDSS 
imaging than in either the ARC or UH88 imaging. While 
the fluxes measured for component A are roughly con- 
stant, component B appears to have faded by ~ 0.8 mag 
between the SDSS and ARC/UH88 epochs. While we 
have already concluded that systematic uncertainties in 

10 For the g and r fluxes we used the V and R values from i]2.3l 
and relaxed the constraint to ±0.25 mag. 
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the photometry may be large (> 0.1 mag), from exper- 
imentation with GALFIT we consistently find that B is 
brighter than A in the SDSS imaging. The relative fluxes 
of the two components may have changed due to intrin- 
sic variability of the source quasar observed with a time 
delay between the two images, and/or differential mi- 
crolensing of the two images. Component B is close to 
the center of the lens galaxy and its observed flux de- 
crease is quite large; therefore, microlensing would seem 
to be the most likely explanation. 

2.5. MMT Spectroscopy 

Spectroscopic observations were obtained on 2009 Nov 
11 with the Red Channel instrument on the MMT 6.5m. 
These observations were nearly contemporaneous with 
the ARC imaging, occuring only five days later. A sin- 
gle 1800s exposure was obtained with a 1.5" slit and the 
270 mm -1 grating centered at 7500A. This provides a 
dispersion of 3.52 A pix at a spectral resolution of 
R ~ 600. The FWHM of unresolved emission lines from 
wavelength calibration exposures is 3.6 pix. The spatial 
pixel scale is 0.297 arcsec pix -1 . The seeing was 1" and 
thin clouds were passing overhead during the observa- 
tion. The instrument was rotated to PA=36.6° in order 
to align the two components from the SDSS imaging in 
the slit. 

Lamp exposures were obtained immediately subse- 
quent to the science observation; however, a standard 
star observation was not obtained until the following 
night, using an identical spectrograph configuration but 
with somewhat better seeing. All images were reduced 
with standard IRAF routines, including a bias correc- 
tion, flat fielding from a continuum lamp, wavelength 
calibration from HeNeAr lamps, and relative flux cali- 
bration from the standard star BD+284211 obtained on 
2009 Nov 12. The science image was corrected for cosmi c 
rays using the L.A. Cosmic routines (|van Dokkumll200ll) . 

The 2D spectra are strongly blended; in particular, 
components B and G are separated by only ~ 2 pixels 
(0'.'6) perpendicular to the dispersion direction (rows in 
the MMT image, see Figure [5]). We adopted a heuris- 
tic approach to extraction of the spectra with all three 
compo nents deblended, similar to that of iPindor et al.l 
(2004). First, the three spectra were treated as Gaus- 
sians fit independently to each row of the image, using a 
modified Levenberg-Marquardt algorithm from the SciPy 
library. We found that allowing all Gaussian components 
to vary and solving for all nine parameters led to clearly 
unacceptable results, even with reasonable initial condi- 
tions. We thus found it necessary to apply strong con- 
straints to the Gaussian parameters in order to produce 
acceptable fits. We assumed that the standard star ob- 
servation provides a model of the spectral trace on a 2D 
image, and then used this model to extract the blended 
science spectra, according to the following process: 

1. A template trace for the 2D spectra was obtained 
from the high S /N observation of the standard star. 
This trace was used to derive the variation of the 
position and width of the spectrum along the image 
rows. The trace positions were fit with a 7th order 
polynomial, but they vary by only ~ 1 pixel across 
the full image. 
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Fig. 5. — Two-dimensional spectra from MMT Red Channel ob- 
servations. The spectral image has been binned by a factor of 8 
along rows for display purposes (the spectral resolution is R ~ 600). 
Approximate wavelength positions along the dispersion axis (rows) 
are labeled. The traces for each component derived from Gaussian 
fitting are shown as dashed lines, while the FWHM of the fits are 
shown at various intervals with a horizontal line. The A spec- 
trum is relatively unblended with the other two, while B and G are 
strongly blended. 

2. Along each row of the science image, a Gaussian 
with three free parameters was fit to the leftmost 
columns of component A (those least affected by 
blending from G). This fit was compared to the 
template trace, and used to derive a translation of 
the template along the image columns to match 
the trace of A, as well as to derive a scale factor 
for the widths corresponding to the difference in 
seeing between science and standard star observa- 
tions. This scale factor combined with the tem- 
plate widths provides an estimate of the widths of 
the point source components A and B. 

3. The fit to A from step #2 was subtracted from 
the image. The same fit was then translated along 
columns to the position of B (as derived from the 
imaging) and subtracted again. The residual image 
is an approximation to the spectrum of G alone. 
Using the template trace again, fitting was per- 
formed on the residual image to derive the width 
of the galaxy spectrum along all rows. 
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Fig. 6. — MMT Red Channel spectra of the components of J0946+1835. The upper panel shows the spectra extracted according to the 
procedure described in 32,51 The spectra have been smoothed to a dispersion of 10. 6A, and scaled to the i— band fluxes derived from the 
image fitting of each component (Table [3J|. Constant offsets have been added to the spectra of A and G for display purposes. The striking 
similarity of the spectra for A and B, including the line profiles and absorption blueward of N V and C IV, provides strong evidence for 
the lensing hypothesis. Prominent elliptical galaxy absorption features arc marked above the G spectrum, using the redshift for the galaxy 
derived from the Ca absorption seen in the SDSS spectrum (Figure[2]l. The A and B spectra have S/N ~ 30 pix -1 on the red side of Lya 
and S/N ~ 10 pix -1 on the blue side. The lower panel shows the ratio of the B and A spectra after smoothing, with no scaling of the 
extracted spectra applied. The locations of the prominent quasar emission lines are indicated. The median ratio is ~ 1.1 in the continuum 
regions, but significantly lower around the emission lines. Uncertainties of ±1(7 in the ratio are shaded around the median value. Note that 
this is the observed flux ratio, and may differ from the actual ratio of the two components due to relative slit losses. 



4. The previous steps provide estimates for the posi- 
tions and widths of each spectral component based 
on expected positions derived from the imaging ob- 
servations. However, the slit was likely not aligned 
perfectly on the center of each component. We then 
used the derived fits as initial conditions to a fit on 
the full 2D image. For this fit, the centers of the 
traces for each of the three components were offset 
by random amounts of up to ±1 pixel along the 
columns. Fitting across the full image was then 
repeated 500 times with different initial offsets, al- 
lowing only the trace positions to vary. In other 
words, copies of the template trace derived from 
the standard star spectrum were varied in position 
until the best alignment with the three components 
was found. The positions from the best \ 2 among 
the repeated fits were then adopted as the best-fit 
positions for the 2D spectra. 

5. Finally, with all positions and widths fixed, a fit 
was performed allowing only the amplitudes to vary 
along each row. The integrated Gaussians for the 
resulting amplitudes and widths were taken to rep- 
resent the ID spectra for each of the components. 

Figure [6] shows the extracted spectra for the three com- 
ponents. The spectra have been scaled such that the syn- 
thetic z-band fluxes derived from the spectra match the 
i-band fluxes from the ARC imaging. 

Several features of the spectra are worth noting. First, 
both the emission and absorption features of A and B 
show remarkable agreement, providing clear evidence 
that they represent lensed images of a single quasar at 



a redshift z = 4.8. This cannot be due to incorrect ex- 
traction of the 2D spectra, as the Gaussians used to ex- 
tract these two components are cleanly separated. While 
the line profiles of A and B are almost indistinguishable, 
the ratio of the two spectra shows that the emission line 
to continuum ratios are different. Component B is well 
within the effective radius of the galaxy, and differential 
microlensing of the continuum a nd broad line regions o f 
B could explain this effect (e.g., Richar ds et al.ll2004af l. 
Alternatively, spectral variability could produce the dis- 
crepant line to continuum ratios, and cannot be ruled out 
by our single-epoch spectra. While the relative fluxes of 
the extracted spectra for the two components are uncer- 
tain due to possible differential slit losses, and the extrac- 
tion was hampered by the strong blending of the B and 
G spectra, we are confident that the observed difference 
between the line to continuum ratios is not an artifact 
of our data. We have experimented with a large num- 
ber of fitting parameters for extracting the 2D spectra, 
and in no viable fit does the B spectrum appear to be 
simply scaled from the A spectrum; in other words, the 
line ratios of the two components are always found to be 
significantly different from the continuum ratios. 

In Section 12.41 we noted color differences between the 
two components of A(i — z) « 0.25 mag, and considered 
whether differential reddening was affecting the observed 
colors. If one of the images, namely component B, was 
strongly affected by extinction from the lens galaxy, this 
would appear in the spectra as a difference in the shapes 
of the continua measured for the two components. In- 
stead, the ratio of the continua measured for the two 
components is flat across ~ 2000A of the observed wave- 
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TABLE 5 
Lens Modeling Results 



Re(") 1.34 ±0.03 
e 0.41 ±0.04 
e (°) -49.0 ±1.4 
At(days) 142.7 ± 6.8 
A'tot 3.2 



length range. However, some color differences should 
be introduced by the observed emission line differences. 
We derive synthetic i- and z-band magnitudes from the 
MMT spectra and find that the emission line differences 
contribute ~ 0.13 mag to the redder color of component 
B. This reduces the net color difference between the two 
components from the ARC imaging to ~ 0.12 mag, which 
is within the uncertainties. 

The absorption features in the galaxy spectrum are 
exactly as expected for a luminous elliptical galaxy at 
z = 0.388, as derived from the Ca absorption features in 
the SDSS spectrum. This provides further confidence 
that the extraction was successful, as neither do the 
galaxy absorption features appear in the quasar spec- 
tra, nor do the quasar emission features appear in the 
galaxy spectrum. 

We determine the quasar redshift to be z = 4.8 from 
the peak of the relatively narrow Lya line, using the 
spectrum obtained for component A. The C IV line has 
significant absorption of its blue wing. A single Gaussian 
fit to the C IV line with the absorbed region masked re- 
sults in a lower redshift of z — 4.77; however, the broad 
component of the C IV line is of ten blueshifted with re - 
spect to the systemic redshift (Richard s et al.l [2002b). 
We fit a power law continuum to the A spectrum and 
three Gaussians to the C IV line, one for a broad com- 
ponent, one for a narrow component, and one for the 
strong absorption. From this fit, the redshift of the nar- 
row component is z = 4.80, in good agreement with the 
Lya redshift, while the broad component is blueshifted 
by - 2000 kms" 1 . 

While the quasar spectra do show absorption features, 
particularly between Lya and N V and blueward of C IV, 
the width of the C IV absorption is not strong enough 
to meet the criteria gener ally adopted for Broad Absorp- 
tion Line (B AL) quasars (jWevmann et al.lll.9911 ) , though 
it has a non- zero Absorption Index (AI = 140 km s _1 ) 
as defined bv lHall et al.l (|2002f ) and may be a mini-BAL. 
This type of absorption appears in a minority of quasars 
and its identical shape in the spectra of both quasar im- 
ages is strong evidence that they are lensed images of a 
single source as opposed to a binary quasar. 

3. LENS MODELING 

We model the system as a Singular Isothermal El- 
lipsoid (SIE) lens generating two i mages of the source 
quasar. We use the lensmodel code {Kccton 20 Oil) to cal- 
culate the properties of both the lensing galaxy and the 
quasar. The positions and fluxes of each of the quasar 
images and the position of the lens galaxy obtained from 
the GALFIT fitting of the ARC z-band images are used 
to provide a total of 8 constraints. We solve for the posi- 
tion and flux of the source quasar, as well as the Einstein 
radius (Re), ellipticity (e, defined as in Table [3]), posi- 
tion angle (8 e ), and position of the lensing galaxy. As 



the model has zero degrees of freedom, the lensmodel so- 
lution converges to x 2 ~ 0. The results derived from the 
ARC i-band imaging are given in Table [5l along with la 
uncertainties (A^ 2 = 1). When fitting the lens model, 
positional uncertainties were increased by 50% and the 
flux uncertainties for all components were set to 5% in 
order to account for systematic uncertainties in deblend- 
ing. The known redshift of the lensing galaxy (z = 0.388) 
enables us to estimate a time delay of 143 days. 

There are some discrepancies between the galaxy light 
profile as measured from the imaging (Table [2]) and the 
derived mass model (Table [5J . If these discrepancies are 
real, one possibility is that they are due to shear (e.g., 
iKeeton et aHll998f) from the nearby galaxies that were 
not included in the model. 

4. DISCUSSION 

Quasars in the SDSS have received extensive follow- 
up to identify lensed objects. The SQLS identifies can- 
didate lensed quasars from the spectroscopically con- 
firmed quas a r sam ple using selection criteria outlined in 
lOguri et al.l ((2006). First, known quasars that are re- 
solved in the imaging may represent two or more lensed 
images at small separations that are blended into sin- 
gle photometric objects within the SDSS database (mor- 
phological selection). Second, large separation lenses are 
selected by searching for nearby objects with similar col- 
ors as the quasa r (colo r selection). Using these meth- 
ods, llnada et al.1 (|2009f ) recently identified five z > 2.2 
lensed quasars in the SDSS. Two are z > 3 quasars 
targeted by the hi gh redshift qua sar selection criteria 
used by the SDSS (jRichards et~aTll2002al) ; for both ob- 
jects, the two quasar images are separated by > 1.5" and 
were successfully deblended by the photometric pipeline. 
J0946+1835 bears a s trong resemblance to the z = 3.6 
lensed quasar found bv lJohnston et ail (|2003l ). which was 
also a galaxy targeted as an LRG by the SDSS and shown 
to be a high redshift lensed quasar by the presence of 
strong, broad emission lines in the SDSS spectrum. 

It is instructive to consider whether J0946±1835 could 
have been selected as a lensed quasar candidate in the 
SQLS even if the lensing galaxy had not been tar- 
geted for spectroscopy by the SDSS. The SDSS registers 
J0946+1835 as two photometric objects, one a point- 
like object representing component A, and the other a 
resolved object representing a blend of B and G. The ob- 
ject representing component A is classified as stellar and 
its colors are relatively unaffected by the lensing galaxy, 
and thus it could have been considered as a quasar target 
by the SDSS. In fact, the object is well within the color 
selection criteria used by the SDSS to target z > 4.5 
quasars (jRichards et al.ll2002af ). with a dereddened PSF 
magnitude of i = 19.52 and colors r — i = 1.82 and 
i — z = 0.06. It was not considered as a quasar target be- 
cause it has the photometric flag DEBLEND_NOPEAK 
set in multiple bands, indicating that no peak was found 
in the deblending. This means that the deblending, and 
thus the colors, are questionable, and the SDSS quasar 
target selection algorithm excludes such objects from its 
color selection. 

Had component A not been flagged DE- 
BLEND.NOPEAK, it would have been a primary 
quasar target in the SDSS. Once identified as a quasar, 
it would be further considered as a lens candidate 
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by the SQLS. Component A has a stellar classifi- 
cation in the SDSS and its PSF likelihood is high 
enough (star_L(r)=0.06 and star_L(i)=0.07) to fail to 
meet the morphol ogical selection criteria of the SQLS 
(|Qguri et al.l I200 6f) . In addition, the strong blending of 
the quasar and galaxy flux in the B+G object causes 
significant color differences between the two SDSS 
objects ((r — i)a — (r — i)b+g = 0.93), such that it 
would also fail to satisfy the color selection criteria. 
Thus, in the hypothetical (though generally expected) 
case that the stellar SDSS object was targeted as a 
quasar instead of the galaxy object being targeted as an 
LRG, this object would not have been identified as a 
lens candidate. It was fortuitous that the lensing galaxy 
was the object targeted for spectroscopy, as the blended 
quasar and galaxy features in the spectrum made the 
nature of the object evident. Furthermore, since the 
B+G object was spectroscopically identified as a quasar 
but clearly has the morphology of a galaxy, J0946+1835 
will be included in the SQLS DR7 additional sample 
by modified mor phological selection (see Section 5 of 
llnada et al.ll2008l ). 

Many of the difficulties encountered with selecting this 
object as a gravitational lens candidate arise from the 
bright lensing galaxy. The galaxy effectively hides the 
second quasar image, both in terms of detecting a sec- 
ond PSF image, as well as the blending of the galaxy 
and quasar colors. The galaxy also affected the well- 
separated quasar image by inducing additional photo- 
metric errors. High-redshift quasars have extremely red 
colors and are detected in only a few optical bands, and 
it is easy to imagine that similar lensed quasars may be 
hiding behind bright lensing galaxies. One way to find 
such objects is in the spectra of galaxy targets, as was 
the case here. This method has been used explicitly to 
find lensed background galaxies in the spe ctra of LRGs 
by th e Sloan Lens ACS Survey (SLACS; iBolton et all 
1200(1 . 

5. CONCLUSIONS 

We have identified the highest known redshift lensed 
quasar. J0946+1835 consists of a lens galaxy at z = 
0.388 and a pair of quasar images produced from a source 
at z — 4.8. An object representing a blend of the lensing 
galaxy and one of the quasar images was targeted as a lu- 
minous red galaxy by the SDSS and classified as a quasar 
based on the spectrum. Subsequent observations show 
that the quasar images are separated by 3" and have 
a relatively low total magnification of 3.2. Comparison 
of photometry from three epochs indicates that there is 
variability of the two quasar images. Deblended spectra 
of the three components show that the two images repre- 
sent a single quasar, and suggest that microlensing may 
affect the relative line equivalent widths, though spectro- 
scopic observations spanning multiple epochs (preferably 
with higher spatial resolution) are needed to confirm this 
interpretation. 

Continued photometric and spectroscopic monitoring 
of this object can be used to measure a time delay to com- 
pare with the val ue derived he re, which constrains the 
Hubble constant (lRefsdallll964D . Multiple spectroscopic 
epochs can also be used to look for variation of the quasar 
emission line equivalent widths consistent with a mi- 
crolensing signature fe.g.. iRichards et aLll2004al) . which 



constrains the size of the broad line region. Lensing can 
also distort the host galaxy image away fro m the bright 
activ e nucleus such that it can be resolved (jPeng et al.l 
12006ft . Probing the host galaxy of this object with HST 
imaging has the potential to allow the relationship be- 
tween black hole mass and bulge mass to be examined 
for the most distant lensed quasar currently known. 

This is the third high-redshift lensed quasar found in 
the spectrum of an object targeted by the SDSS as an 
LRG, joining the z = 3.6 lensed quasar reported by 
IJohnston etall fl2003 |) and a z = 2.2 len s from the SQLS 
high-redshift sample (|Inada et al.ll2009h . The SDSS-IlO 
will target 1.5 million LRGs to a flux limit ~ 1.5 mag 
fainter than the SDSS-I/II; this fifteen-fold increase in 
the number of LRG spectra has the potential to result in 
a similar increase in discoveries of high-redshift quasars 
lensed by those galaxies and revealed by their spectra. 

Another appro ach to identify qua s ars is by their 
varia bility (e.g., ISchmidt et al.l I2010t IKozlowski et al.l 
120101). Planned syn optic surve ys such as PAN- S TARRS 
([Kaiser et aHmm and LSST (jlvezic et al.ll2008l) will be 
able to identify high redshift quasar lenses through sys- 
tematic searches for pairs of stellar objects at small sepa- 
rations that both vary (with a time dela y), or by looking 
for re s olved objects (galax i es) th at vary (jKochanek et al.l 
I2006D . IQguri fc Marshalll (|2010fl predict that thousands 
of new lensed quasars will be discovered by upcoming 
time-domain surveys, including several hundred at z > 4 
from the LSST alone. 
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